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Abstract 

We consider a multidimensional cosmological model with FRW type metric having 
4-dimensional space-time and (i-dimensional Ricci-flat internal space sectors with a 
higher dimensional cosmological constant. We study the classical cosmology in com- 
mutative and GUP cases and obtain the corresponding exact solutions for negative 
and positive cosmological constants. It is shown that for negative cosmological con- 
stant, the commutative and GUP cases result in finite size universes with smaller 
size and longer ages, and larger size and shorter age, respectively. For positive cos- 
mological constant, the commutative and GUP cases result in infinite size universes 
having late time accelerating behavior in good agreement with current observations. 
The accelerating phase starts in the GUP case sooner than the commutative case. In 
both commutative and GUP cases, and for both negative and positive cosmological 
constants, the internal space is stabilized to the sub-Planck size, at least within the 
present age of the universe. Then, we study the quantum cosmology by deriving the 
Wheeler-DeWitt equation, and obtain the exact solutions in the commutative case 
and the perturbative solutions in GUP case, to first order in the GUP small param- 
eter, for both negative and positive cosmological constants. It is shown that good 
correspondence exists between the classical and quantum solutions. 

PACS numbers: 98.80.Hw; 04.50. +h 



1 Introduction 

It is well known that quantum description of gravity must be considered when we want to 
deal with systems in the Planck scales, such as very early universe or strong gravitational 
field of a black hole. In the absence of gravity, quantum description of a system can be de- 
rived from classical consideration by replacing of Poisson bracket with usual commutation 
relations ({ , } — >■ 4[ , ]). When we want to consider the gravitational effect in quantum 
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description of a system, some essential modification in the ordinary quantum principal are 
needed. General Uncertainty Principal (GUP) is a modification of Heisenberg Uncertainty 
Principal in the Planck scale. Such a generalization has already been considered in the 
context of string theory where the string can not probe distance smaller than the string 
size PQ-jH]- Some general view to the GUP and its application to cosmology are proposed 
in [7J- [13] and jH], [15], respectively. Michele Maggiore in the discussion of a Gedanken 
experiment for the measurement of the area of apparent horizon of a black hole in quantum 
gravity , using rather general and model independen consideration, showed that a mini- 
mum length of order of the Planck length in a Generalized Uncertainty Principal emerges 
naturally from any quantum theory of gravity [TO]. In [7J, the simplest form of the GUP in 
a one dimensional system, is written as: 



where Lpi = y — 10 is the Planck length and j3 is a positive constant to be of 
order unity. The new term in the above equation is important when x,Ax ~ Lpi. It is 
possible to show that the above GUP relation (equation [1]) can be derived from the following 
generalized Heisenberg algebra: 



Since cosmology provides the ground for testing physics at high energy, it seems natural 
to expect the effects of quantum gravity in this context. Alternatively, in cosmological 
systems, the scale factor, matter field and their conjugate momenta play the role of dy- 
namical variables of the system; so, introducing GUP in the corresponding phase space is 
particularly relevant. 

In the past few years the search for a consistent quantum theory of gravity and the quest 
for a unification of gravity with other forces have both led to a renewed interest in theories 
with extra spatial dimensions [H]-[T9]. Extra spatial dimensions ideas date back to the 
seminal work of Kaluza [20] and Klein [21] . These extra spacial dimensions must be hidden 
and assumed to be unseen because they are compact and have small radius, presumably 
with typical dimensions of the order of planck length, O(10~ 33 cm). At Planck time, t p = 



= O(10 44 s), the characteristic size of both internal and external dimensions are likely 



to have been the same, and the internal dimensions may have had a more direct role in the 
dynamics of evolution of the Universe. 

As we discussed above, the existence of extra dimensions and the influence of GUP 
become evident at high energy. So, it is clear that in studying cosmology at high energy, 
we should presumably consider extra dimensions and GUP together [22], [23] . 

The paper is organized as follows. In section 2, we construct a multidimensional cos- 
mology. By convenient coordinate transformation, we calculate it's Hamiltonian which 
describes an isotropic oscillator-ghost-oscillator system. In section 3, we investigate the 
equations of motion in classical model, at first in commutative case and then by consid- 
ering the influence of GUP. Section 4, devotes to quantum cosmology. We construct the 
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Wheeler-DeWitt (WD) equation in commutative and GUP cases and solve it to obtain the 
wavefunction of the corresponding universes. 



2 The Cosmological Model 

It is well known that universe in large scale is homogenous and isotropic and has 3- 
dimensional space at least in the order of experimental tests. So it is possible to exist 
unseen compact internal space (Extra dimension) with small radius. For our investigations, 
we consider a mult i- dimensional cosmological model in which the space-time is established 
by a FRW type metric of 4-dimensional space time and a d-dimensional Ricci-flat internal 
space [24] : 



ds 2 = -dt 2 + - R2<y !\. (dr 2 + r 2 dQ 2 ) + a 2 {t)gfMx l dx\ (3) 



where k — 1, 0, —1 represents the usual spatial curvature of external space in FRW model, 
R(t) and a(t) are the scale factors of the external and internal space respectively, and 
is the metric of the internal space which is assumed to be Ricci-flat. The total number of 
dimensions is D = 3 + d. The Ricci scalar can be derived from the metric© 
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where a dot represents differentiation with respect to t. We consider an Einstein-Hilbert 
action functional with a D-dimensional cosmological constant A: 

S = WTtI d D x^(n-2A)+S YGH , (5) 

where ko is the D-dimensional gravitational constant and Sygh is the York-Gibbons- 
Hawking boundary term. Substitution Eq.fllJ), after dimensional reduction we have: 

S = -v D „ x J dtl^6R 2 ^R + 6MR 2 + ^-^^R 3 -6k<l>R + 2<l>R 3 A^, (6) 

where 

-(=)■ 

and ao is the compactification scale of the internal space at present time. We can set 
U£>_i = 1. To make the Lagrangian suitable, we consider the following change of variables^- 

&R 3 = T 2 (xl - xl), (8) 



1 Here we have some differences from |24| in denning the parameters of equations (|10j) . 



$ p+R t- = T{ Xl + x 2 ), 
= T(xx-x 2 ). 



(9) 



with 




2 V d + 2 



(10) 



where R = R(xi,X2) and $ = ^(xi,x 2 ) are functions of new variables Xi,x 2 . Using the 
above transformations and concentrating on k = 0, the Lagrangian becomes 

C = (x 1 *-x^) + ±(±^)(xl-xl). (11) 
We can write the effective Hamiltonian as@ 

Pi , , ,2J2\ ( V\ , , ,2„2 



where 



«= + (12) 



Equations (|T2|) . (|T3|) and (jSJ) show that the potential energy for our oscillator-ghost- 
oscillator system [29] is proportional to vacuum energy A times the volume of the multidi- 
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mensional universe and the total momentum of such a system is pf ot = y — 

The effective oscillators in Eq.(lT2l appear to be decoupled and the dynamics of two 
dynamical variables, namely X\ and x 2 seems to be independent of each other. However, 
the Hamiltonian constraint (H = 0) imposing on oscillators connects their dynamics 0. 
Therefore, the dynamics of external large dimensions R and internal compactified extra 
dimensions a are inevitably related to each other by the equation ([9l). 



2 In the Hamiltonian of cosmological models having matter contents, usually a wrong sign between 
geometric and matter sectors appears. This is due to the fact that the matter content of the universe 
has positive energy while the gravitational content has negative energy. In the present model, there is no 
explicit matter content, however, we may interpret the sector of extra compact dimensions as a fictious 
matter field. Then, the Hamiltonian of the variables R, a having a sign difference between the four and 
extra dimensional sectors is transformed to that of the variables xi , x-i with the same sign difference between 
X\ and X2 parts. These Hamiltonians are usualy called oscillator-ghost-oscillator. 

3 We know that general relativity is a time reparametrization invariant theory. Every theory which is 
diffcomorphism invariant casts into the constraint systems. Therefore, general relativity is a constraint 
system whose constraint is the zero energy condition {% =0). 



3 Classical solutions 



3.1 Commutative case 

As in [SUES], the dynamical variable defined in ([9]) and their conjugate momenta satisfy: 

{Xp,Pv}p = 11l<v, ( 14 ) 

where rj^ u is the two dimensional Minkowski metric and { , }p represents the Poisson 
bracket. The equations of motion can be written as: 

PtM = {p^U} P = -2lj 2 x^ (15) 

Combining two equations in (Tl5|) we obtain 

x fJ , + u) 2 x li = 0. (16) 

First, we assume a negative cosmological constant. According to (1131) . u> 2 is then pos- 
itive. Considering f fl2|) . it is obvious that Eq. (ll6j) describes the equations of two ordinary 
uncoupled harmonic oscillators whose solutions are 

x,(t) = A^ t + B^ t , (17) 

where and are constants of integration. Imposing the Hamiltonian constraint 
("H = 0) introduce the following relation on these constants. 

A^B" = 0. (18) 

Finally, using (J7]) and (Q, the scale factors take the following forms 

a(t) = k\ [sin(wt + 0i)] d ^+"+-p-"-) [ s in(u;t + 4> 2 )} d{p+ a +- p - a - ) 



R(t) = k 2 [sin(ujt + (f) 1 )} p +' T +- p -' T -[sm(ut + (f) 2 )] p +°+~ p - a - , (19) 

where k\ and k 2 are arbitrary constants and (pi and (p 2 are arbitrary phases. If we consider 
the Hamiltonian constraint, the following relation is imposed on these constants 



d + 3 



kfk% cos(0i - 2 ) = 0. (20) 



Because of ki, k 2 ^ 0, this leads to (pi — (p 2 = | . In what follows, we investigate the behavior 
of a universe with one internal dimension (D = 3 + 1). 
By setting (pi = f and (f) 2 = 0, we find: 



R(t) = k 2 a/ sin(a;t), 

a(t) = h (21) 



whose classical loci is given by 



a = 1V1 - i? 4 , (22) 

where the values of k\ and &2 are absorbed in the values of R and a. Figure. 2 (left), 
shows this loci for negative cosmological constant. It is seen that the dynamics of R is 
constrained by the dynamics of a and vice versa. This is expectable, because the dynamics 
of two variables is subject to the zero energy condition so called Hamiltonian constraint. 
Using these solutions, we can calculate the Hubble and deceleration parameter for both 
R(t) and a(t) as 

H R {t) = ^ = r ot(ut), 

= _ R(t)R(t) =1 + 2tan 2 (a;t) 

H a (t) = = ~(cot(ut) + 2t&n(ut)), 

a(t) 2 

a(t)d(t) 2cos 2 (u;t)(5 + cos(2u;t)) 

Qa[t) ~ (-3 + cos(2u;t)) 2 ' ( } 

From ()2ip . we see that as time increases from t = toward j-, R(t) and a(t) are 
increasing and decreasing functions of t, toward the maximum and zero, respectively. At 
t — j-, R(t) reaches to its maximum and a(t) is contracted toward zero. From t — j- 
to t — -, R(t) is contracted toward zero while a(t) is decreasing in the region of negative 
values. Although the negative value for a(t) is meaningless, but it is a 2 (t) > which 
is physically viable in the definition of metric0. Therefore, we stick to the behavior of 
a 2 (t) > instead of a(t). As is seen in Fig.l (solid lines), from t = j- to t = -, both 
R 2 (t) and a 2 {t) become decreasing and increasing functions, respectively. According to the 
behavior mentioned above, the four dimensional sector begins from a Big Bang at t ~ 0, 
then expands till t — j- toward a maximum value after which the gravity takes over the 
expansion and this sector starts contracting toward a big crunch at t ~ -. In order to 
avoid of initial singularity, we will assume the Big Bang to occur at the Planck time tpi. 
In fact, it is reasonable that at planck time we set R(tpi) = a(tpi). In this regard, the 
extra dimensional sector begins from a 2 (tpi), then contracts till t — toward a vanishing 
minimum value after which this sector starts expanding toward a 2 (^ — £pz) = a 2 (tpi). we 
see that during the whole time evolution of universe (tpi < t < - — tpi), the scale factor 
of internal space is contracted toward zero and can never exceed a(tpi). By considering 
present value of Hubble constant, we see that the age of universe is t present = -cot _1 (^ a ) 
whose order of magnitude is in agreement with present observations, namely u^ 1 ~ 10 17 s. 
So, the present universe is in tideway to get to the maximum and minimum of R 2 {t) and 
a 2 (t), respectively, within At ps 0.57a; -1 . 

In the time interval - < t < — , R 2 (t) and a 2 (t) become negative Figs.l (solid lines). 
To avoid this problem we may consider two options: 



4 This situation is similar to that of wavefunction in quantum mechanics, where it is the square of the 
wavefunction which is physical and not the wavefunction itself. 
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Figure 1: Time evolution of the (squared) scale factors of universe with one extra dimension 
and negative cosmological constant. Solid and dashed lines refer to the scale factors in 
commutative and GUP framework respectively. Left and right figures are the external and 
internal dimensions respectively. 



1) Classically, negative values for R 2 (t) and a 2 (t) are nonphysical, so we may leave 
the cosmology with imaginary scale factors as the nonphysical solutions of the Einstein 
equations and think that the physical universe will end at t — -, with no further extension 
or history. 

2) Quantum mechanically, as the scale factor R(t) goes to zero (when t approaches 
the big crunch singularity), the strong quantum gravitational effects can drastically change 
the initial conditions so thati? 2 (t) and a 2 (t) may become again nonzero positive functions 
capable of establishing another Big Bang. 

By, considering 

R(tpi) = k 2 ^/sm(utpi), 

a(t Pl ) = h \ Pl) (24) 

the initial condition R(tpi) = a(tpi) results in 

^ = 10 61 , (25) 
h 



by which we obtain the following ratio 

= 10 ei tan(wt). (26) 



"(0 , M M 



a(t) 

If we consider the present radius of external space to be equal to the radius of observed 
universe, namely 10 28 cm, then we see that the present radius of internal space is about the 
Planck length (10 _33 cm). 

By considering smallness of the internal dimension and using (T251) to evaluate the ap- 
proximate magnitude of H a (t present ) ~ uj ~ 10~ 17 , we see that the variation of a{t) at 
present time is ignorable. So, in the context of negative cosmological constant, the com- 
pactification and the stabilization of the internal space at present status of the universe is 
established. 

Now, we assume a positive cosmological constant for which u 2 is negative. If we replace 
u 2 with — u 2 in Eq.( fl6l) . the new solutions are then obtained by replacing trigonometric 



0.0 0.5 



1.0 1.5 2.0 2.5 3.0 3.5 



Figure 2: Left figure is classical loci in the negative cosmological context and right figure 
is the classical loci in the positive cosmological context. 

functions by their hyperbolic counterparts in the solutions (j!9p . By considering Hamiltonian 
constrain we obtain the following solutions 



a(t) = 
R(t) = 

where for d = 1 we have 



ki [cosh(cjt)] d (p+°+-p-°-) [sinh(wt)] d W CT +-p- CT -> 

-p- p+ 
k 2 [cosh(ujt)} p +' T +- p - 



a(t) 



[sinh(wt)] p +' T +~ p - 
cosh(cot) 



(27) 



y/sinh(u;£) 



whose classical loci is given by 



R(t) = k 2 \/ sinh(wt), 



(2? 



(29) 



Figure. 2 (right), shows this loci for positive cosmological constant. As in the negative 
cosmological constant case, the dynamics of both R and a is constrained by the Hamiltonian 
constraint. The scale factor ratio is obtained 



R{t) 

a(t) 



10 bi tanh(u;t), 



(30) 



and we have 



H R {t) 

q R {t) 
H a (t) 

Qa(t) 



R(t) co . . 



_R(t)R(t) 
a(t) co 



1 - 2tanh 2 (wt), 



a(t) 2 
a(t)a(t) 



— coth(ojt) + 2tanh(u;£)), 

2 cosh 2 (cot) (5 + cosh(2u;t)) 
(-3 + cosh(2wt)) 2 ' 



(31) 
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Figure 3: Time evolution of the scale factors of universe with one extra dimension and 
positive cosmological constant. Solid lines refer to the scale factors in commutative case, 
dashed and dot-dashed lines refer to the scale factors in GUP framework with /3 = 10 -4 
and (3 = 10~ 5 respectively. Left and right figures are the external and internal dimensions 
respectively. 

Note that, like the case of negative cosmological constant, the magnitude of the radius of 
external to internal space is asymptotically (t — > oo) about 10 61 . Equations (|28|) show that 
R(t) is an increasing function of time. But, a(t) at first decrease with time till t ~ 0.88a; -1 
and then increase exponentially (see Fig. 3). As we mentioned above, at planck time, the 
characteristic size of both internal and external dimensions are assumed to be the same, 
namely (R(tpi) = a(tpi)). If we consider the age of universe about u^ 1 ~ 10 17 s, we see 
that at present time we are around the minimum point of a(t). Also by some calculation 
on a(t), we find that in the time interval tpi < t < lAluo" 1 , a(t) can never exceed it's initial 
value at Planck time a(tpi). Therefore, the internal scale factor remains small for a very 
long period with a duration of about 140 times of the present age of the universe. 
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Figure 4: Left and right figures are respectively Hubble an deceleration parameters of 
external space for universe with one extra dimension and positive cosmo logical constant. 
Solid lines refer to the scale factors in commutative case, dashed and dot-dashed lines refer 
to the scale factors in GUP framework with (3 = 10 -4 and = 10~ 5 respectively. 

3.2 GUP case 

In more than one dimension, it can be shown that the generalized Heisenberg algebra 
corresponding to GUP is defined by the following commutation relations [7J [5] 



[x h pj\ = i(S i:j + /35 i:j p 2 + (3'piPj), (32) 

where p 2 = YliViPi an d /3(3' > are considered as small quantities of first order. Throughout 
the whole process we work in the units with H — G — c — 1. We assume that momenta 
commute with momenta 

\pi,Pj] = 0. (33) 

Using Jacobi identity [[xi, Xj],pk] + [[xj,pk], Xi) + [[pk,Xi],Xj] = the commutation relation 
for the coordinates are obtained as: 

fa, x i\ = i— — Tf^^ ~ ~ ^ vix 3 - Pj x *y ( 34 ) 

It is clear from the above equation that the coordinates do not commute. So, we can 
not work in the position space to construct Hilbert space representations. But, if we choose 
the special case j5' = 2(3, we can see from equation (|34|) that coordinates commute to first 
order in (3 and thus we can work in coordinate representation. By the following definitions, 
equations ( |32l) and ( |33l with f3' = 2(3 can be realized up to first order in j3 (see Appendix) [30J 



Xi = x i0 , pi = p i0 (l + /3pl), (35) 
any Hamiltonian can be written as [30] 



where [xio,pjo] = iSij, p\ = ^PioPio and pio = —ig^- We can show that the p 2 term in 
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Figure 5: Left and right figures are respectively hubble an deceleration parameters of inter- 
nal space for universe with one extra dimension and positive cosmo logical constant. Solid 
lines refer to the scale factors in commutative case, dashed and dot-dashed lines refer to 
the scale factors in GUP framework with (3 = 10 -4 and (3 = 10 -5 respectively. 



p 2 = v\ + 2(3p 4 . 



(36) 



In our model, by considering /3' = 2/3 and in first order in /3, the commutation relation 
between position and momentum operators can be summarized as 



[x uPl ] = i(l + (3p 2 + 2(3pl), [x 2 ,p 2 ] = i(l + Pp 2 + 2(3p 2 2 ), 



(37) 



[xi,P2] = [x 3 ,pi] = 2ifip x p 2 , 



(38) 



} = [Pi,Pj}=0, i,j 



1,2. 



(39) 



where p 2 tot = \(jp\—p 2 )- We aim to investigate the effects of the classical version of GUP. As 
we know, we must replace the quantum mechanical commutators with the classical poisson 
bracket as [P,Q] —> i{P,Q}. Thus, in classical phase space the GUP deformed poisson 
algebra is achieved from (!3"7l) - (!39l by considering [P,Q] — >■ i{P,Q}. 



{x 1>Pl } = (1 + (3p 2 + 2f3pj), {x 2 ,p 2 } = (1 + (3p 2 + 2{3p 2 2 ), 



(40) 



{x 1 ,p 2 } = {x 2 ,pi} = 2fip x p 2 , 



(41) 



{xi,Xj} = {pi,Pj} = 0,i,j = 1,2, 



(42) 



In [31] such a deformations algebra is used on classical orbit of particles in a central force 
field and on kepler third law. It is notable that this modification is significant at Plank scale 
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Figure 6: Left and right figures show the classical trajectories in %\ — x 2 plane for a universe 
with one extra dimension in the context of negative and positive cosmological constant, 
respectively. 

and we need quantum description. But, before quantizing the model we want to construct a 
deformed classical cosmology. Note that, in transition from quantum commutation relations 
to Poisson brackets we Shall keep the GUP's parameter j3 fixed as h — > 0. The equations 
of motion can be written as (see |32j ) 

x x = { Xl ,H} = ^ Pl (l + 5Pp 2 ), (43) 



x 2 = { X2 ,H} = --p 2 (l-3f3p 2 ), (44) 

Pl = { Pl ,H} = -2co 2 x 1 (l + (3p 2 + 2f3p 2 1 )+4uj 2 ^x 2 p 1 p 2 , (45) 

P2 = {P2, H} = 2u 2 x 2 (l + (3p 2 + 2(5p 2 2 ) - Atu 2 (3x lPl p 2 . (46) 

These show that the deformed classical cosmology forms a system of nonlinear coupled 
differential equations that can not be solved analytically. 

In the context of negative cosmological constant, u 2 is positive ( fl3|) . Numerical solutions 
of Eqs. (|43l - P6l) shows that R(t) and a(t) like the classical ones have periodic behavior, but 
the time interval between Big Bang and big crunch is shortened (dashed lines, Fig.l). For 
later period, the time interval becomes longer with respect to first one while the maximum 
value of R(t) is shorter. Comparison of the results in the present GUP framework with the 
previous classical commutative case (Fig. 1) reveals that in the present framework we have 
larger external and smaller internal dimensions. Also the rate of variation of a(t) is less 
than the previous case. So, the compactification and the stabilization of internal space in 
GUP framework is more favored than the classical commutative case. 

By replacing oj 2 with — oj 2 in Eqs. (j43]) - (l46|) . we can get the corresponding equations in 
the case of positive cosmological constant. By numerical analysis, we find that at early 
time the scale factors behave like classical commutative case (Fig. 3), but at late time both 
the external and internal dimensions are greater than those of commutative case. Also, 
investigation in time dependence of the Hubble and deceleration parameters shows that at 



early time these parameters nearly behave like classical commutative case, but at late time 
they behave very differently. 

Some recent observations indicate that the universe is currently undergoing an acceler- 
ating period of expansion. This is consistent with the results obtained here considering a 
positive cosmological constant. Depiction of the Hubble and deceleration parameteres for 
the scale factors R and a can be of useful help to understand the origin of current accel- 
eration in this multidimensional model. To this end, we have depicted Hr, and H a , q a 
within the figures 4 and 5 (see solid lines). Figure 4 shows that becomes negative a little 
bit earlier than ut ~ 1, namely the present age of the universe. This means, the accel- 
eration of the universe has started recently in the present multidimensional commutative 
case. On the other hand, Fig.5 shows that q a is always negative and has a minimum at the 
position where q^ becomes negative. Therefore, it seems the behavior of q a is responsible 
for the behavior of q^ and vice versa, as is interpreted in the following. The present model 
describes a multidimensional cosmology, with a positive cosmological constant. Typically, 
a standard 4-dimensional FRW cosmology with a positive cosmological constant predicts 
an accelerating behavior of the scale factor. In a multidimensional cosmology, however, 
it is reasonable to think that the overall repulsive force due to the positive cosmological 
constant manifests as an interplay between the accelerating and decelerating behaviors of 
the external and internal scale factors. Looking at the figures 4 and 5 reveals that, typ- 
ically for both commutative and GUP cases, at the beginning of time, qn is positive (R 
is decelerating) while q a is negative (a is accelerating). As time passes, qn approaches the 
threshold of negative values (R is less decelerating ) while q a goes to more negative values 
(a is highly accelerating). When q^ enters the region of negative values (R is accelerat- 
ing) the q a reaches its minimum (a stops the increasing acceleration; the minimum is not 
shown in Fig.5). Finally, q^ becomes more negative (R is highly accelerating) and q a goes 
to rather less negative values (a is slowly accelerating). In order to better compare the 
results obtained in commutative and GUP cases, we have depicted all the GUP diagrams 
for two values of the parameter /3 with a difference of ten times order of magnitude. It is 
seen that as the value of ft becomes smaller, the results of GUP more coincide with the 
results of the commutative case. Considering a very small value for the parameter (3, for 
example suggested by string theory, reveals that at least at the present age of the universe 
(t ~ w _1 ) it is impossible to distinguish between the GUP and commutative cases in that 
which approach fits better the data. 

4 Quantum solutions 

4.1 Commutative case 

To investigate the quantum version of the model, we use the WD equation, "H^ = 0. Here 
"H is the operator form of the Hamiltonian (Tl2|) . Using the canonical procedure to quantum 
mechanics that the phase space variables are replaced with quantum operators, by replacing 
Pi = ~ ^lix * n we get the following WD equation: 




(47) 



This equation is the quantum version of oscillator-ghost oscillator system with zero energy 
condition. We can solve it by Formal variable separation approach as: 



^ni ,ri2 

{x 1 ,x 2 ) = U ni {xi)V n2 {x 2 ), (48) 



d 2 W- 

+ (A - Au 2 xl)Wi = 0, Wi{i = l,2) = U,V (49) 



in negative cosmological background, the solutions are 



\ V 4 2 

U ni {xi) =[ — ) -==H ni (V2^x 1 ), (50) 

7T I y/2 ni ni\ 



1/4 



K 2 (x 2 ) = — -===^(^2^), (51) 
7T / 2 n. 2 ! 

with A = 2rii + 1 and the restriction n\ — n 2 — n. Here, H n (x) are Hermite polynomial 
and the eigenfunction are normalized as: 

/+oo 
e~ x2 H n (x)H m (x)dx = 2 n n 1 / 2 n\5 nm . (52) 
-oo 

The general solution of the WD equation can be written as a superposition of the above 
eigenfunctions: 

* 2 ) = (— ) V2 e-^ + ^ V ^ 1 H n (V2^x 1 )H n (V2^x 2 ). (53) 

n 

Coefficients c„ are so chosen as to make the states coherent 1 33 1 : 



i 

4 ,nn\ 



- ^-e-i'-l 2 . (54) 
2u) (_i)T n j 

In general, the quantum solutions do not offer semiclassical description of some space- 
time region unless we introduce a decoherence mechanism which is widely regarded as 
necessary to assign a probability for the occurrence of a classical space-time. However, 
in the lack of desired decoherence mechanism, in order for a satisfactory classical-quantum 
correspondence is achieved we may at least investigate if the absolute values of the solutions 
of Wheeler-DeWitt equation have maxima in the vicinity of the classical loci. Fig. 7 (left) 
shows the square of the wavefunction in the case of negative cosmological constant. As we 
can see from this figure and the corresponding classical trajectory in Fig. 6 (left), the peak 
follows the classical trajectory. Thus, it is seen that there is a good correspondence between 
quantum description and the classical pattern of the model. 



Figure 7: The square of the wavefunction of the universe in x\ — x^ plane with one 
extra dimension and negative cosmo logical constant. Left and right figures refer to the 
commutative and GUP framework, respectively. 



We get the corresponding answers in the positive cosmological background by replacing 
uj 2 with —uj 2 in the above equations. Fig. 7, shows the square of the wavefunction. Because 
of highly increasing dependence of the amplitude of the wavefunction on x\ and X2, the 
fluctuation of wavefunction on the line x\ = £2 can not be seen using the normal scales 
adopted in this figure. However, the peak of squared wavefunction follows properly the 
classical trajectory, namely the line X\ = x^ in Fig. 6 (right). Therefore, similar to the 
case of negative cosmological constant, there is a good correspondence between quantum 
description and the classical pattern. 



4.2 GUP case 

Here, we aim to investigate influence of GUP in quantum cosmological model that has 
been presented above. The Hamiltonian of the model is given by (Tl2|) . To construct WD 
equation in GUP framework we use Eq. (l36l) with the representation pn) = —i-Q§-^ for the 
momentum operator. In first order in /3, by ignoring zero subscript we have 
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*(ari,x 2 ) = 0. (55) 



We see that this is not a separable differential equation, so we set a procedure to 
convert it into a separable one. Up to first order in /3, we assume a perturbative solution, 
$ = ^r(o) + /^M, w here # (0) is the wave function of commutative case, namely the solution 
of Eq. pTl) . Substitution of this perturbative solution into Eq. fl55|) reveals that due to the 
presence of j3 in the first order part of Eq. fl55|) we may put ^ ~ hence we can use the 
commutative case equation ( |47|) for this part. By this consideration we have 
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*(xi,x 2 ) = 0, 



(56) 
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*(a?i,ar 2 ) = 0. 



(57) 



Because of the term x\-ks, the equation is inseparable too. To resolve the problem we 



use the classical loci x\ + x\ 



(see figure 6). 



pA pA pp Q2 \ f pp pp. 

£ ° 8u 2 (r 2 - x 2 )— + 8a; 2 — -(x 2 ) + ( -— + — 

dxf dx\ " 2 dx\ dx\ 2 J \ dxf dx\ 



¥ = 0. 



Separation of the variables as ^{x\, X %) = U{x\)V{xz) with perturbative solutions in 
first order of (3 as U{x l ) = U^{ Xl ) + pU^{xi) and F(x 2 ) = V {0 \x 2 ) + (3V^(x 2 ), and 
perturbative separation of the separation constant A = A^ + (3\^ leads to 



in zeroth order in j3 



d2W ^ + (A (°) - 4a, 2 * 2 )iy(°)(^) = 0, (5f 
ox i 

W^( Xi ;i^ 1,2) ^U^(x 1 ),V ( -°\x 2 ), 



in first order in /?. 

+ (A (0) - Auj 2 x 2 )W^\ Xi ) = g( Xi ), (59) 
WW{ Xi ;i = 1,2) = uW( Xl ),V<»(x 2 ), 

We can see that Eq.( |58|) is commutative limit (/? — > 0) of Eq. (l55|) (see Eq. (|49l ). So it's 
solution is 



\ V 4 2 

Wffo) = — -==ff n (^). (60) 
7T / V2™n! 

Eq.( )59l) is a inhomogeneous differential equation. The solution of homogeneous part is 




Figure 8: The square of the wavefunction of universe in x\ — x 2 plane with one extra 
dimension and positive cosmological constant. 



= Wi°\ Xl ). (61) 
By use of Eq.(16Tl). the solution of inhomogeneous part can be written as 



(I) = U^)Wi 0) (x t )+^)Wi°\x t ), (62) 



where 



1 f Xi g n Wn^ 

Uxi) = ^w^wf-wrw^ 1 

ct s r ^nW^%-g n W^ A 

KAXl) J wWwF-wWwW *' 

. f n 2nW$»W<!» . \ .... 

and a prime denotes for differentiation with respect to x$. Combination of zeroth and first 
order solutions leads to 



W n ( Xi ) = W<?\ Xi ) + P(W^( Xl ) + U^W^(x t ) + ^(xi)W^( Xi )). (64) 

Therefore, we may write the general solution of WD equation in the GUP framework as a 
superposition of above eigenfunctions 

^gup(x u x 2 ) = ^2c n W n (x 1 )W n (x 2 ). (65) 



We set, as in Eq. (154j) . the coefficient c n to the coefficient of commutative wavefunction. 
Figure. 7 (right) shows the square of the wavefunction with negative cosmological constant 



in GUP framework. As is seen in this figure, there are four peaks distributed around 
Xi,X2 = ±8.5 which in the classical viewpoint correspond to a state with R max . Also, 
the overall peak follows the classical trajectory in figure. 6 (left). Thus, similar to the 
commutative case, there is a good correspondence between quantum description and the 
classical pattern of the model in GUP case. 

5 Problem of time 

The canonical approach at the classical level gives the Hamiltonian constraint. Promoting 
this constraint to the quantum level does not give a time-dependent wave equation, rather 
leads to a stationary and timeless equation so called Wheeler-DeWitt equation. This sug- 
gests that noting at all evolves in the universe. Therefore, one is faced with the problem 
of explaining the origin of the notions of time in the laws of physics. This is so called time 
problem in quantum gravity and there are different approaches to solve this problem [26J. 
Among these approaches, a new and relevant approach to this paper is the one introduced 
in [27]- [28] where the authors offer an explanation of time evolution in quantum gravity 
by using of the probability density in quantum mechanics. They establish a method called 
"Probabilistic Evolutionary Process" (PEP) which states that if the system is in a given 
quantum state, then a small perturbation inspires this state to end up in another quantum 
state with the higher probability of existence. The evolutionary process of initial quantum 
states toward the more probable final quantum states suggests a possible solution to the 
problem of time in quantum cosmology. 

In our problem by consideration of time evolution in classical context, which accords 
with increasing in x\ and X2, we may say that such a transition in quantum state from 
low probability to high probability, can be considered as a quantum evolutionary process. 
This interpretation of time evolution may be justified in the positive cosmological constant 
case, where there is considerable difference in the probability amplitude of wavefunction. 
However, because of nearly equal probability amplitude of wavefunction in the negative 
cosmological constant case, this interpretation of time evolution is not of practical use to 
solve the time problem. 

In classical context, the time evolution is in accordance with going from R m i n to R ma x- 
As we mentioned in previous subsection, by consideration of higher probability of existence 
in quantum mechanics and using the discussion in references [27J-[28], we may explain 
the time problem of the WD equation in quantum gravity. It should be noted that, this 
interpretation is due to the influence of GUP. 

6 Conclusion 

We have considered a multidimensional cosmology having FRW type metric of 4-dimensional 
space-time and (i-dimensional Ricci-flat internal space coupled with a higher dimensional 
cosmological constant. The classical cosmology in commutative and GUP cases are stud- 
ied in detail and the corresponding exact solutions for negative and positive cosmological 
constants are obtained. For negative cosmological constant, both cases result in finite 
size universes which differ by their size and age, while, for positive cosmological constant 



both cases result in infinite size universes having late time accelerating behavior in good 
agreement with the current observations. Both commutative and GUP cases with negative 
and positive cosmological constants result in the stabilization of internal space to the sub- 
Planck size. We have also derived the Wheeler-DeWitt equation and obtained the solutions 
in both cases for both negative and positive cosmological constants. It is shown that good 
correspondence exists between the classical and quantum solutions. 

It is worth noting that the effects of GUP in this model are important not only in the 
early universe but also in the late time behavior of the cosmic evolution. So, in the GUP 
framework, we can see somehow the indirect quantum gravitational effects at large scale. 
Moreover, the influence of GUP to solve the time problem in this model is remarkable. 

7 Appendix 

We want to show that the momentum pj in the GUP framework in terms of the ordinary 
low energy momentum pj , equation ( )35|) fulfill the commutation relations ( 1321) and ( |33l) 
with (3' = 2(3 up to first order in /3. Let us compute ( 132|) 



where we have used the identity [P, QS] = Q[P, S] + [P, Q]S. Up to first order in j3 it is 
obvious that 



[a>j , Pj\ 



[xio,Pj + PPjoPkoPko] = [xi ,Pjo\ + (3[x i0 ,p j0 pkoPko} 
iSij + (3pjo[xio,PkoPko\ + (3[x i0 ,p j0 \pkoPko 
i5ij + 2i/3S ik pkoPjo + iPSijPwPko 
i(Sij + f35ijpl + 2f3p i0 p j0 ), 



(66) 



Pp 2 



/3piPi = f3pio + f3p i0 PkoPko = fiPioPio = PpI 



(67) 



and similarly /3piPj = (3piPj. Substituting these in equation (1661) yields 



[xi,pj] = i(8ij + (36ijp 2 + 2(3piPj), (68) 
which is nothing but Eq. (l32l) with f3' = 2(3. Eq. ( |33l) is a direct consequence of [pio,Pjo] = 0. 
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